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Abstract 

The (ordinary) Sachs-Wolfe effect relates primordial matter perturbations 
to the temperature variations 5T/T in the cosmic microwave background ra- 
diation; 5T/T can be observed in all directions around us. A standard but 
idealised model of this effect leads to an infinite set of moment-like equations: 
the integral of P{k) j'^{ky) with respect to A: (0 < < oo) is equal to a given 
constant, Ci, for £ = 0, 1, 2, . . .. Here, P is the power spectrum of the pri- 
mordial density variations, je is a spherical Bessel function and y is a positive 
constant. It is shown how to solve these equations exactly for P{k). The same 
solution can be recovered, in principle, if the first m equations are discarded. 
Comparisons with classical moment problems (where jj{ky) is replaced by 
k^) are made. 
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I. INTRODUCTION 



Some might say that the uhimate inverse problem is to understand the origin of structure 
in the Universe using information that is currently available: it is the central problem in 



early-Universe cosmology ■ One of the available cosmological observables is the cosmic 

microwave background radiation (CMBR). It is known that the temperature of the CMBR is 
remarkably uniform in its spatial variation, with 5T /T ~ 10~^-10~^ [^|, |18, §1.5], implying 



that the expansion of the Universe is largely isotropic. Conversely, inhomogeneities in the 
density of the Universe lead to temperature anisotropics, so that these can be used as a 
sensitive test of theories of structure formation IRB . 



In these theories, the primary unknown function is -P(fc), the power spectrum of the 
primordial density fluctuations. Under certain assumptions, it can be shown that P{k) 
satisfies the following set of equations: 

/ k-^P{k)jj{ky)dk = C,, £ = 0,1,2,.... (1) 

Here, y is a given positive constant, ji is a spherical Bessel function and the constants Ci 
are given. 

Experiments have provided estimates for a finite number of the C^'s. How can these 
be used to recover P{k)7 This is a major question, but it is not our main concern here. 
We are interested, first, in an idealised problem: given an exact knowledge of all the C/s, 
reconstruct P. We show that this inverse problem can be solved exactly. This is a new 
result. Furthermore, we show that the same solution can be recovered if the first few (in 
fact, any finite number) of the Ces are unknown. 

What use is an exact inversion formula for an idealised problem? First, it reveals the 
ill-posed nature of the problem.^ Thus, P{k) is found to be given by the Fourier sine 
transform of a certain function go{^), which is defined by an infinite series with a finite 
interval of convergence. This means that techniques of analytic continuation are required, 
a process that can be very difficult numerically [|T2| and one that gives the problem its ill- 
posed character. Second, our inversion formula reveals some of the analytic structure. For 
example, the low-A; behaviour of P{k) is intimately related to the asymptotic behaviour of 
5fo(A) as A ^ 00. Third, exact results can be used to test numerical algorithms designed for 
the finite-data problem. 

The plan of the paper is as follows. In the next section, we sketch a derivation of the 
governing equations (|I|). Careful derivations can be found in the literature; references are 
given. Our aim is to motivate the study of (|I]) in a way that is accessible to non-specialists. 
Thus, we limit our discussion to perhaps the simplest model of the underlying physics. 
In section we formulate two moment-like problems, called the Basic Problem and the 
Reduced-data Problem in which we are given Ce for i > and i > io, respectively, where io 
is any fixed positive integer. Apart from the system (^, we also consider (in section |V|) the 
analogous classical moment problem where jj{ky) is replaced by k^. For both cases, we apply 



^For the definition of an ill-posed problem, see, for example |8|, ch 4]. 
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a general method described in section This replaces the infinite system of moment-like 
equations by a single integral equation. Application to the astrophysical problem is made 
in sections |VIHIX|. 



II. THE GOVERNING EQUATIONS 

There are several physical causes of temperature variations in the CMBR (Kolb and 
Turner list five ||T8| , p 383]). For large angular scales, the dominant contribution comes 
from variations in the gravitational potential; this is known as the Sachs-Wolfe effect 
An elementary discussion of the relevant equations can be found in [l^, §9.6.2], p5| , §6.4] 
and in p8H7|,|6| , whereas advanced, detailed, derivations that are up-to-date with the present 
state of computation can be found in |jl5| and ||2^. A brief derivation is sketched here; for 



a rigorous derivation, involving time- dependent processes within General Relativity, see the 
cited books and papers. 

Consider the standard case of a flat universe {Qq = 1) with zero cosmological constant 
(A = 0). The observed temperature fluctuations of the CMBR in the direction of the unit 
vector y are 

STjy) _ T{y)-To 
T To ' 

where Tq is the mean temperature. The Sachs-Wolfe effect relates the observed ST/T to 
the primordial density fluctuations in the early Universe at a time when light and matter 
decoupled. For the microwave radiation, the whole Universe at this time is referred to as 
the last scattering surface 0, p 74]. Since the time of last scattering to the present, the 
distance travelled by light is y = 2c/ Hq, where c is the speed of light and Hq is the Hubble 
constant. A present-day observer at a given point will receive light from a spherical shell of 
present-day radius y on the last scattering surface. Thus, it is natural to expand on that 
shell in terms of spherical harmonics y^(^), 

^-^ = Y.^Uy)Yriy): (2) 

£,Tn 

where aim are dimensionless expansion coefficients. 
We can also give a Fourier representation for 6T/T, 



T 



V J 0(fc) exp{iyk-y)dk, (3) 



where Q{k) is a Fourier amplitude. The factor V = V/{2tt)^ where V is a (large) volume, 
so that G is dimensionless. (For simplicity, we have used a continuum description (Fourier 
transform) here rather than a discrete description (Fourier series) followed by an appropriate 
limit.) 

We can compare (|) and (|^), using the plane- wave expansion 



exp {iyk ■y)=A7rJ2z^ ^ky) Y^iy) Yr{k), 
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where je is a spherical Bessel function, k = kk, k = \k\ and the overbar denotes complex 
conjugation. We find that 



aUy) = 4vr/V J &{k) je{ky)Yr{k) dk, 
whence use of the addition theorem for Legendre polynomials gives 

^ = 47r(2£ +l)V' ff Q{k) Q(l7)j,{ky) ^k'y) P,{k • k) dk dk' . (4) 

The Sachs-Wolfe effect relates 5T /T, through physics on the last scattering surface, to 
the primordial density fluctuations at time of last scattering, 5p/ p. The latter quantity has 
a Fourier representation similar to namely 



P 



V J A{k) exp{ik-x)dk. 



It is usual to assume that A (A;) is an independent Gaussian random variable of zero mean. 
Thus (A(fc)) = and 



[A{k) A(fc')) = V-^P{k) 6\k - k'), 
where P{k) is the power spectrum of the density fluctuations. Note that 

P{k) = V J (A(fc) A(fc')) dk'. (5) 

Then, the Sachs- Wolfe analysis shows that (the Fourier transforms of) ST/T and Sp/p are 
related by 

A{k) = Ak^Q{k), (6) 

where ^ is a real constant with the dimensions of area. (Roughly speaking, the temperature 
fluctuations are proportional to the perturbed gravitational potential, $ say, and is 



proportional to the density fluctuations; from ||28|, we have A = —2c^/Hq = —^y"^-) So, if 
we take an ensemble average of (^ and deflne rotationally-invariant dimensionless coefficients 

Cehy 



1 



2i 



|2 



tm 



we find that 



Ci = ^J k-^P{k)j^{ky)dk 

POD 

= {4n/AfV k-^P{k)3j{ky)dk. (7) 

JO 
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The set of coefficients {Cg} is known as the angular power spectrum.^ 

The temperature fluctuations have been measured experimentally; for example, this was 



done by the COBE satellite (see |^5[, and references therein). This gives an estimate of 
the angular power spectrum 0, §7]: thus, we have values of for a range of £- values. From 
these, we would like to extract the power spectrum, P{k). This is a moment-like problem. 
In section [Vll| , we give an explicit solution of an idealised problem, in which we know Ce for 
all i>0. 

In practice, there are complications. First, the measured values of Cq and Ci, the 
'nuisance parameters' 0, p 175], are unreliable. We show that P{k) can be recovered, 
nevertheless, using Ce for all i > io, where Iq is any fixed positive integer. 

Another complication is that we cannot actually measure Ce because it is defined as an 
ensemble average: we always have a relative uncertainty (or cosmic variance) of (2£+l)~^/^, 
even if the experimental measurements are perfect [@], p 354], p 160], @, p 189]. 
Once inverted, these errors will lead to corresponding errors in the calculated P{k). As the 
errors in Ce decay as i increases, one possibility is to use Ce for i > io, where £q is chosen 
sufficiently large, so as to recover P. One could then compute Ce for i < Eq, and compare 
with the experimental estimates, giving the possibility of cross-validation. However, we defer 
discussions of computational aspects to a later paper. 

Further complications come from the underlying physics. Thus, the Sachs-Wolfe effect 
typically is decomposed into two additive quantities, the 'ordinary' and the 'integrated' 
Sachs- Wolfe effects, of which (|^) only represents the former. The methods developed in this 
paper are only applicable to the ordinary Sachs- Wolfe effect. In the idealized limit that the 
CMBR anisotropy is composed only of super-Hubble scale, adiabatic, density perturbations 
in a fiat {Qq = 1) universe, the ordinary Sachs- Wolfe effect dominates. This ideal limit 
has an approximate validity in inflationary cosmology, which is a promising candidate for 
development as a fundamental theory of the early universe; for reviews, see, for example, p!8| , 
ch 8] or p5| , ch 10]. The understanding of the CMBR anisotropy is still a developing area 
of research, although the general consensus is that there are substantial modifications to 
the idealized regime of the ordinary Sachs- Wolfe effect. We will not attempt to survey 
this vast body of research here, but the interested reader can begin with a few exemplary 
papers |1l6| , |T0| , p^ , |l^ |7| . It is safe to say that, when applied to real data, there would be 
practical limitations to the methods developed in this paper. Nevertheless, they offer a new 
and very different perspective on the subject 



^In General Relativity, a density fluctuation is not a gauge-invariant quantity. In practical terms 
this means density fluctuations are not well-defined observables. For sub-Hubble scales, the gauge 
dependence of density fluctuations is negligible, so that any reasonable gauge choice is operationally 
acceptable. However for super-Hubble scales, the gauge dependence is acute. Properly, the power 
spectrum P{k) in (^) must be defined with respect to 'gauge-invariant' quantities ||5[^ that in the 
sub-Hubble scale regime are readily identified as density fluctuations. This is a well-studied problem 
in the literature [24,1£ ,26,32, p^, 2 1,29, p!4|. For our purposes, however, these details are unimportant 
and we will be content with our slightly imprecise description, dH). What is important for our work 
is that the basic relation between P{k) and the angular power spectrum Ce is of the form (|^). 
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III. FORMULATION OF TWO MATHEMATICAL PROBLEMS 



Let us begin by considering the following moment-like problem. 

Basic Problem. Given a set of numbers {C„} and a set of functions {ipnik)}, with 
n = 0, 1, 2, . . find a function f{k) (in a certain space) that satisfies 

/ fik)Mk)dk = Cn, for n = 0,1,2,.... (8) 



We are mainly interested in one particular set {ipn}, given by 

Mk)=fn{ky), (9) 
where y is fixed. For a simpler model problem, we shall also consider the choice 

Mk) = k\ (10) 



which gives rise to a classical (Stieltjes) moment problem |33,^. 



We are also interested in the following related problem where the first m of are 
omitted. 

Reduced-data Problem. Given a set of numbers {C„} and a set of functions {tpnik)}, 
find a function fm{k) that satisfies 

fmik) Mk) dk = Cn, forn>m>0. (11) 

Here, m is fixed and fo = f. 

It is important to note that we seek in the same space as /. We are also using the 
same numbers {C„} and the same functions {ipn} (for n > m) as in the Basic Problem: 
thus, Cn and ipn are independent of m. 

As we saw in section ||, the Reduced-data Problem, with tpn{k) = j^{ky), arises in 
connection with the Sachs- Wolfe effect. In that context, the basic unknown function is the 
power spectrum P{k) = k'^f{k), apart from a constant factor; see (^. 

IV. A GENERAL METHOD 

A general method for treating moment-like problems is to replace them by an integral 



equation. Thus, choose a second set of functions {0„(A)}, multiply the n-th equation of (|TT]) 
by (f)n and sum over n. This gives 



/•oo 

/ fm{k)lCm{k,\)dk = gm{X), (12) 
Jo 

where 

oo 

/C^(fc,A)= E^nW<^"(A) (13) 
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and 



9m{X) = E Cn0n(A). (14) 

n=m 

To be effective, one has to be able to solve the integral equation (|12[), perhaps by recognising 
the left-hand side as a known integral transform; implicitly, this requires that the sum 
defining the kernel Km can be evaluated in closed form. To make the method rigorous, one 
has to show that the sums in (0) and (|l^ are convergent, and that the implied interchange 
of summation and integration is justified. Note that, at this stage, A is unspecified. 

We will use this method below for the two sets {V'n} given by ([T0|) and (^. We remark 
that a similar method was used in p2| , pB| to analyse the so-called null-field equations, a 
system of moment-like equations (involving Hankel functions) that can be used to solve the 
boundary-value problems for acoustic scattering by bounded obstacles. 



V. MOMENT PROBLEMS 

We start with ipn{k) = k^, giving the following simple Reduced-data Problem. 
Moment Problem Aim- Find a function fm{k) that satisfies 

roo 

/ fm{k)k''dk = Cn for alln>m>0. (15) 
Jo 

We assume that fm is a smooth function that has a Laplace transform. 
We can rewrite ( [T3| ) as 

POO 

/ {k"'fm{k)}k^dk = Cn+m, forn = 0, 1, 2, . . .. 
Jo 

Multiply by </)„(A) = {—X)^/n\ and sum over n to give 

roo 

/ {k'^fm{k)}e-^'dk = gm{X), (16) 
Jo 

where 

oo / \\n 

9m{\) = E Cn+m^^. (17) 

n=0 

Equation (|16]) gives the Laplace transform of k"^fm{k), so that we can invert to obtain fm 
itself. This shows that fm is obtainable uniquely from the given coefficients {Cn} with n > m. 

In the derivation above, we assumed that the series (|17[) is convergent. In fact, we 
only require convergence for |A| < Aq, say, and then define gmW for larger |A| by analj^ic 
continuation. The assumption that fm{k) has a Laplace transform is enough to ensure that 
k"^fm{k) also has a Laplace transform, this being what is actually required in the derivation. 

The connection between classical moment problems and the Laplace transform is well 
known; see, for example, p 97] and 0, p 230]. For a connection with a Hankel transform 
(take 0„(A) = (-A)V(n!)2), see |3|, p 96]. 
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How are fm and fo = f related, for m > 0? From (0), we have 



(n-m)! \ dA; „t;; " n\ 



m 



We also have 



I {k"'fm{k)}e-'''dk=i-—\ f,n{k)e-^Uk. 
So, integrating m times gives 

roo 

/ Uk)e-^Uk = go{X)+Prn{X), (18) 

JO 

where po = and 

m—l 

Pm(A) = ttiX'^ 

e=o 

is an arbitrary polynomial of degree (m — l). But the left-hand side of ([ISD is a Laplace 
transform, and so it must vanish as |A| ^ cxd (in a right-hand half-plane). go{X) has the 
same property, as it is the Laplace transform of /q. Hence, the polynomial PmW must be 
absent: 

roo poo 

/ U{k)e-^''dk = go{X)= h{k)e-^Hk. 
Jo Jo 

Thus, fm{k) = fo{k) for all m. This means that if we know, a priori, that fm and /o are 
both in the same space (here, we assumed that they both have Laplace transforms), then 
they are equal: deleting the first m moments C„ does not represent a loss of information. 
Similar comparisons can be made between and fm' with m ^ m' . 
For a simple example, take C„ = n\ for n > 0, whence 

oo 

^?o(A) = EMr = (1 + ^)"' 

n=0 

for |A| < 1. Note that we can define go{\) for all A 7^ —1 by analytic continuation. By 
inspection, we have f{k) = e~^. 



VI. THE SACHS WOLFE EFFECT 

Here, we assume that ipnik) = jl{ky), where jn is a spherical Bessel function and y is 
fixed. Thus, we consider the following problem. 

Sachs— Wolfe Problem Vm- Find a function fm{k) that satisfies 
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/ fmik) jliky) dk = for alln>m>0. (19) 
Jo 



Here, m is fixed, and the constants Cn are given. 

What conditions on are appropriate? Let us assume that 



k^~^ as /c — > oo. 



fmik) 

Then, considering the integral on the left-hand side of (|I9|) , we see that we need 

2m + a + 1 > for convergence at /c = 

and 

P > for convergence at k = oo. 

Note that convergence at = depends on the smallest value of n used, which is m. 
However, we want conditions that do not depend on m. So, given that the Basic Problem 
is Vo, minimal conditions are 

a > -1 and ^ > 0. (20) 

Arguments based on causality imply that P{k) = 0{k'^) as k ^ P|,PD|], where P{k) = 
k'^ f{k). Thus, the desired physical solution should have a >2. This has led to studies of 
model power spectra of the form 

^^^^ = l + (^/fco)-+/3-i' ^^^^ 
where a > 2, /3 > and /cq is a constant [|]. 



Next, let us consider an exact solution [|ri|, formula 6.574.2]: 

-2.. , _ 7r/-^r(/j)r(n- 1/3 + 1) 

^ ^"^^^^ ~ 2/^+Hr(i(/3 + i))pr(n + i/5 + i)- ^^^^ 

This holds for < (3 < 2n + 2. Note that the right-hand side is 0{n~^) as n — * oo. 

This exact solution throws some light on the next question: how does the integral on 
the left-hand side of (|l^) behave as n — >■ cx)? If we split the range of integration at k = 1, 
say, we can easily see that the asymptotic contribution from integrating over < A; < 1 is 
exponentially small. The dominant contribution comes from integrating over > 1, so that 
the large-fc behaviour of fm{k) can be used. Comparing with the exact solution given above 
shows that 



Cn = 0{n ^) as n 



oo. 



Thus, there is a link between the large-n behaviour of the (given) coefficients C„ and the 
large- behaviour of fm{k). 

Another exact solution is [O, formula 6.577.1]: 
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i"00 ^ 
Jq ^2 ^ ^2 " 2^ In+l/2{hy) ^n+1/2 ( ^^ol/) , (23) 

where /„+i/2 and -ft'„+i/2 are modified Bessel functions and is a positive constant. Note 
that the right-hand side is 0{n~^) asn ^ oo. This formula gives the angular power spectrum 
for the model (^Tj), exactly, when a = 2 and (3 = 1. Further formulae can be obtained by 
differentiation of p3|) with respect to k^. 

VII. THE BASIC PROBLEM 



We are going to apply the general method of section |I^ to problem Vq, which we restate 
here. 

Problem Vq. Find a function f{k) that satisfies 

POO 

I f{k)jl{ky)dk = Cn for n = 0,1,2,.... (24) 

We have to choose the set {<^„(A)}. There are very few known sums involving products 
of Bessel functions. From Abramowitz & Stegun formula 10.1.45], we have 

N „ / 9/, N sin \ky\/2 — 2 cos 6} 

^ (2n + 1) P„(cos 6) jliky) = } (25) 
n=o kyV2-2cose 

where P„ is a Legendre polynomial and 6 is unrestricted. If we define A by 

A = yV2 - 2 cos ^ = 2y sin ^6, (26) 

we can rewrite (^5]) as 

,^ , , sin kX 

K:oik,\) = ^^, (27) 

where /Cm is defined by (|1^) and 

0„(A) = (2n + 1) P„ (l - liX/yf) . (28) 
Note that, for convergence of the integral in (|12|), we require that A is real. Moreover, we 



can obtain the formula (27) for all A > if we allow 6 = Or + iOi to be complex; specifically, 



we can suppose that 9 lies on the L-shaped contour given by 

{0 < < TT with 6*7 = 0} and {0 < 6*/ < cx) with Or = tt}. 
So, multiplying (^^ by (2n + 1) P„(cos0) and summing over n gives 

roo 

/ A;"V(A;) sinAA;dA; = A^o(A), for A > 0, (29) 
Jo 

where 
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goiX) = E (2^ + 1) Cn Pn (l - \{\lyf) ■ (30) 



n=0 

This reduces the determination of f{k) to the inversion of a Fourier sine transform: 

IT Jo 



f[k) = — A^o(A) sinAfcdA. (31) 

vr Jo 



The derivation assumes that k~^f{k) has a Fourier sine transform, which is consistent with 
the convergence conditions (|20|). It also assumes that the series (^) is convergent for some 
values of A, and that go{X) can be defined for other values of A by analytic continuation. 

For an example, take C„ to be defined by the right-hand side of (|2^) for n > with 
< /5 < 2. In order to calculate go{X), defined by (pOD, we use the following expansion (see 
the appendix for a derivation): 



{1 - x)-^ = Y.{2n + 1) Cn Pn{x) for |x| < 1, (32) 

n=0 



where 7 < 1 and 



~7 



r(l-7) r(n + 7) 



(33) 



r(7) r(n-7 + 2)- 

Comparing with C„, we set 7 = 1 — |/3 whence 

Cn = {A/n)2-^/Y'^r{(3) sin (ivr/^) C„. 
With X = 1- liX/yf, we find that 

gJX) = ^- (34) 

' ^ ' 2 r(/3) sin (ivr/?) ^ ' 

Note that the series defining g^^X) converges for < A < 2?/, but we can use (|^) to define 
go{\) for all A > 0. Hence (0) gives 

k 

f{k) = —J- ^ / A'^-^ sin Xk dX = k^-^, 



T{p) sin (ivr/?) Jo 

in agreement with 



VIII. THE REDUCED-DATA PROBLEM Vm 

For m > 0, we consider the Reduced-data Problem Vm (the Sachs- Wolfe Problem) 
described in section |VJ. 

Proceeding as for Vq, we multiply (|1^) by (2?2 + 1) P„(cos^) and sum over n from n = m. 
This gives, using (pU]), 
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fUk) EV"Win(%)| dk = gM), (35) 



where and 0„(A) are defined by ( p^ and (^), respectively. 

Next, multiply ( pS]) by A and write as 

/■oo m— 1 

/ fc-Vn^(A;)sinAfcdfc-A5]0„(A) / /„(A;)j2(%)dfc = A(7™(A). (36) 

Noting that P„(x) is a polynomial in x of degree n, we see that the second term on the 
left-hand side of (^) is a polynomial in A of degree (2m — 1). So, applying the differential 
operator d^'^/dA^™ gives 

^2m »oo (^2m ^2m 

k'^ fm{k) sinXkdk = —— (A^m(A)) = — - (A^o(A)) . 



Integrating 2m times then gives 

/■oo 

/ k-'fmik)smXkdk = XgoiX)+P2mW, (37) 



where p2m is an arbitrary polynomial of degree (2m — 1). But the left-hand side of (^) is the 
Fourier sine transform of k~^ fm{k). This is assumed to exist and, by the Riemann-Lebesgue 
lemma, it must vanish as |A| ^ cxd. A5'o(A) has the same property, as it is the Fourier sine 
transform of k~^ f, by (pOf). Hence, the polynomial P2m(A) must be absent. It follows that 
fm{k) = f{k) for all m. 

The solution given above is correct but somewhat deceptive, for in the formulation of Vm 
we are not given C„ for < n < m. Thus, we cannot form the series (|3^) defining 5'o(A). 
So, instead of fBTI), we obtain 



/ k'"" fmik) sinXkdk = Xgm{X) +p2mW, 
Jo 



where the polynomial p2n{X) is to be determined by the requirement that the right-hand 
side vanishes as |A| — > oo. Thus (just as for Vo), we have to effect the analytic continuation 
of gm{X). Numerical consequences remain to be explored. 



IX. DISCUSSION 

We have given an explicit formula for solving the Basic Problem Vq. For the application 
we have in mind, all the constants Cn are non-negative: does it follow that the solution 
f{k) is positive? In general, the answer is 'no'. We show this by giving an explicit counter- 
example. Let 

f{k) = y{{ky/nY/'-{2MA}, (38) 

where A is a constant to be chosen. Note that if A > 0, f{k) < for /c < AA^/liry). (The 
various factors in (^8]) are merely inserted for algebraic convenience.) We observe that the / 
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defined by ([38| ) is a linear combination of two of the exact solutions (p2D, corresponding to 
/? = I and P = 1. Thus, 



23/2 {r(|)} Yin 



We have to show that C„ > for all n > 0; we do this using an inductive argument. 
Straightforward calculation shows that 

_ 4n + 3 4{n + l) 

4n + 5 " (4n + 5){2n + l){2n + 3) ' 

So, Cn > for all n > provided A> and Co > 0. But 

Co = 1 - A, 

so we can take any A with < A < 1. 

APPENDIX 

Here, we derive the expansion (|32|) with (^). From (^), orthogonality gives 



c„ = i / ^(1 - x) ^Pn{x) dx 



2 



— Tl" /■! d" 

(l-x)-^-j^(l-x2)"dx 



2"+%! 7-1 ' da;*" 
using Rodigues' formula. Integrating by parts n times, using 



[i-xy = -^—^{i-xy 



dx'^' ' r(7) 

and noting that the integrated terms vanish, gives 



2'^+%! r(7) 7-1 
In the integral, put x = 1 — 2t; it becomes 



1 r(n + 7) /\l_3;2)n^^_^^-,-n^^_ 



2n+l-7 1^ ^-7(1 _ = 2"+^" 

7o 



-7 



n!r(l -7) 



/o ' ' r(n-7 + 2)' 

using the integral definition of the Beta function, and the result ( p3D follows. 
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